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A charged Higgs pair can be produced at an ee collider through a t-channel exchange of a heavy 
neutrino (N) via e + e~ — > H + H~ and, if N is a Majorana particle, also via the lepton number 
violating (LNV) like-sign reaction e ± e ± — > H ± H ± . Assuming no a-priori relation between the 
effective eNH + coupling (£) and light neutrino masses, we show that this interaction vertex can 
give a striking enhancement to these charged Higgs pair production processes. In particular, the 
LNV H~H~ signal can probe N at the ILC in the mass range 100 GeV Z m N ~ 10 4 TeV and 
with the effective mixing angle, £, in the range 10 -4 £ 2 £ 10~ 8 - well within its perturbative 
unitarity bound and the PPov limit. The lepton number conserving (LNC) e + e _ — > H + H~ mode 
can be sensitive to, e.g., an O(10) TeV heavy Majorana neutrino at a 500 GeV International Linear 
Collider (ILC), if £ 2 > 0.001. 



PACS numbers: 14.60.St,13.15.+g,13.66.Hk,12.15. Ji 

The discovery of neutrino oscillations, which indicates 
mixing between massive neutrinos, stands as the first 
direct evidence for physics beyond the Standard Model 
(SM). Thus, the SM, or for that matter also the Minimal 
Supersymmetric Standard Model (MSSM) and any other 
multi-Higgs doublet model, has to be expanded to include 
massive neutrinos that mix. Since there is still no under- 
standing of the nature of these massive neutrinos, i.e., 
Majorana or Dirac-like, the extensions to these models 
can go either way. However, from the theoretical point 
of view, the more natural and therefore more appealing 
way is to generate sub-eV Majorana neutrinos by adding 
heavy right-handed neutrino singlet fields (N with a mass 
mjv) and relying on the seesaw mechanism which yields 
m v ~ m 2 D /mjs[, where m v is the solar and atmospheric 
neutrino masses and mu are the Dirac mass terms. Thus, 
the classic seesaw mechanism links light neutrino masses 
(m„ ~ 10 -2 eV) with new physics (mjv) at the GUT scale 
(if mo ~ tuw) or at the multi-TeV scale (if tod ~ m e ) - 
both are well motivated theoretically. In either case, this 
minimal seesaw setup yields a vanishingly small heavy- 
to-light neutrino mixing: Uw °c y m^/uiN, in which 
case heavy neutrinos effectively decouple. 

In this work we will assume no a-priori relation be- 
tween the couplings and light neutrino masses. Indeed, 
interesting models have been proposed, which assume 
that the underlying physics for sub-eV neutrino masses 
goes beyond the standard, minimal, seesaw framework 
(see e.g., pi 0,0, HI), thus allowing the mixing angle U^n 
to depart from its classical seesaw value (i.e., y m„/mjv)- 
Our primary purpose in this paper is to propose col- 



lider experiments that can sensitively probe the relevant 
masses and mixing angles. 

The fact that a Majorana mass term violates lep- 
ton number by two units, i.e. AL = ±2, has strik- 
ing phenomenological implications, some of which have 
been extensively studied in the past few decades - the 
most noticeable one being neutrinoless double beta decay 
(A, Z) — ► (A, Z + 2) + e~ + e~ With the upcoming 
CERN Large Hadron Collider (LHC) and the planned 
high-energy ee machines, the study of LNV collider sig- 
nals in processes involving heavy Majorana neutrinos is 
particularly well motivated. Indeed, such high-energy 
LNV signals have been extensively studied in t he p ast 
two decades: at pp and pp collisions [1, 0, 0,0, [3, at 
an e + e~ [rllllj]. e~e ~ [121 and e~7 collisions [rlll3| and 
at an ep maclhne i I 111 . Il4| . In addition, LNV mediated 
by exchanges of a heavy Majorana N can be manifest in 
top-quark and W-boson decays pl| and in rare charged 
meson decays [8| . 

In this letter we propose that charged Higgs pair pro- 
duction at an ee machine like the ILC (International Lin- 
ear Collider) can serve as a very sensitive probe of heavy 
neutrinos (N) with masses up to ms ~ 10 4 TeV. The 
two processes we will explore are: 



(1) 
(2) 



'Electronic address: atwood@iastate.edu 
^Electronic address: shaouly@physics.technion.ac.il 
* Electronic address: soni@bnl.gov 



These processes are generated at the tree level via a t- 
channel (and u-channel in the case of the like-sign H^H^ 
signal) exchange of a heavy N, which can have a sizable 
coupling to e — H + when » m\y (see below). Note 
also that, while the H + H~ signal is not sensitive to the 
nature of N, i.e., be it Majorana or Dirac particle, the 
like-sign H^H^ signal is generated only if N is a Majo- 
rana neutrino. 
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FIG. 1: The LNV and LNC cross-sections a and a~ + in units of oqed = o(e + e~ — > 7* — > fi + n~), as a function of tun 
and with £ — miii(0.1, £u) or £ = min(l, £u), for a 500 GeV (right figure) and 1000 GeV (left figure) ee collider with charged 
Higgs masses mu = 150 and 300 GeV, respectively. 



In any model with more than one Higgs doublet, a 
right handed neutrino (be it Majorana or Dirac-like) can 
couple to a charged Higgs via its Yukawa interactions 
with the scalar doublets YL^N. The H + £N Yukawa- 
like coupling can, therefore, be written as: 



£>H+W — 



( 3 ) 

where m d and mi are the neutrinos and charged leptons 
Dirac mass terms, respectively, and the /„'s are some 
functions of the Higgs doublets VEV's. For example, 
fv — (fv')~ 1 = tan/3 within an extended supersymmet- 
ric or a non-supersymmetric type II two Higgs doublet 
model with heavy right handed neutrinos. 

Since we will be interested in the effects of a heavy 
Majorana-type neutrino, we assume N to have a Majo- 
rana mass term of the form MNN/2. Thus, we rewrite 
the H + eN interaction Lagrangian ^ in a generic form, 
in terms of the physical mass of TV, ttin, and assuming 
f%m e « f^m D : 



eN 
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^N(l 



2\/2 mw 



j 5 )eH + + h.c. 



(4) 



where the dimensionless parameter £ defines an effective 
H + eN coupling strength that will be the key factor which 
controls the size of the cross-section for reactions (Q} and 
(|2"|). In particular, we absorb in £ any relation between 
mo and win which will depend on the mechanism respon- 
sible for the generation of neutrino masses. Depending 
on the size of £, the interaction term in ((4]) can have strik- 
ing implications for N — H + phenomenology at an high- 
energy ee collider, in particular, for the charged Higgs 
pair production mechanisms in ([T]) and 



As mentioned earlier, we will assume no a-priori rela- 
tion between mjy and light neutrino masses and mixing, 
thus, adopting a model independent approach towards 
the effective H + eN coupling £, and treating it as a free 
parameter (assuming for simplicity that it is real). As 
such, £ should be subject only to existing experimental 
constraints and perturbative unitarity if applicable (see 
below). However, since there is no direct experimental 
constraint on £ that we know of, we will use the seesaw- 
like framework as a guide for estimating the experimen- 
tally allowed size of £. In particular, within a seesaw- 
like mechanism for generating light neutrino masses, the 
Dirac mass mo is related to the physical heavy Majorana 
mass by mr) ~ Uin^n, where Uin ~ mo/mN being the 
heavy-to-light neutrino mixing angle (see e.g., Pilaftsis in 
Q). Thus, in this case we obtain (see Eqs. [3] and [3]): 



U. 



eN 



(5) 



It should be emphasized again that, although the min- 
imal seesaw framework restricts the mixing Um to be 
vanishingly small, one cannot rule out Uin ~ O(l) even 
for mpf ~ 0(100) GeV, based on seesaw-like scenarios 
which go beyond the minimal framework, such as radia- 
tive seesaw, double seesaw and fine-tuned relations or 
symmetries between the mass matrices mjj and tojv, as 
was noted for example in |l|, |2|, y, |j| . 

Using the seesaw- like motivated relation in ((5]) , the ex- 
isting experimental constraints on U c n can be applied 
as an estimate for the allowed size of £. In particular, 
U e N is constrained from precision electroweak data (e.g., 
from invisible Z-decays) \U eN \ 2 < 0.012 at 90% CL 
if mjv > mz- This does not rule out £ ~ 0{1), since f„ 
which depends on the details of the scalar sector of the 
underlying theory, can be of order 10. A stronger con- 



3 



straint on U e N comes from neutrinoless double-/? decay 
{PM E3: J2 N U e N / m N < 5 • 10- 5 TcV- 1 . We note 
though, that the /3/?oi> bound, although tighter, is subject 
to cancelations if there are more than one heavy N's that 
couple to an electron, in particular if they have different 
masses and mixings (see e.g., discussion by G. Belanger 
et al. in 0). 

As noted above, another constraint can be applied 
on the £ — toaj- plane from perturbative unitarity by 
demanding that Tm < toat/2. In particular, in 
our simplified framework N can decay via: N — > 
W^t^, Zvi, Hv£, H^t^, where the decay rates into 
the gauge-bosons is cx | C/ e iv | 2 , see e.g., [13]. Taking 
Mfj » m 2 v ,m 2 z ,m 2 H a,m 2 H+1 the total decay width of 
N (i.e., summing over the above decay channels) is: 

r N « 3|c| 2 r^ • (1 + JL) ; ee £-?fm N . (6) 
V \.fv\ J 647rm^ 

Thus, taking \f^\ 2 >> 1 (i.e., assuming that the total N 
width is dominated by its decays to the Higgs-bosons), 
perturbative unitarity imposes the following bound: W 

t < fo „*/f.»!S„ 0.715^. (7, 
g V o m,N raw 

This gives e.g., £ < 1 for m N > 700 GeV or £ < 0.1 for 
to n ~ 7 TeV. In what follows, we will always apply the 
unitarity bound in iff]), or restrict £ to be smaller than 
either 1 or 0.1 when the unitarity bound does not apply. 

The cross-section for the LNV e~e~ — > H~H~ mode 
is given by: 

a a . BP r%( 2 

(8) 

where (tqbd = cr(e + e — > 7* — > ) = 4na 2 /3/s, 

sw = sin6w, (3 = VI ~~ 4rff, w = 1 + 2r^ — 2rjf and 
r x = m 2 /s. Thus, in the limit m 2 N » s, we obtain 
a cx £ 4 m^y /rriyy, i.e., the cross-section grows with mjv 
for a constant effective mixing £. This growth behavior 
is, however, "cured" when we impose the perturbative 
unitarity bound in |7]): 

,2 m 327r/3 1 

cr (mjy » s)|^ P T . (9) 

In Fig. [T] we plot the LNV cross-section cr~~ in units 
°f vqed, as a function of ton for both a 500 and 1000 



[1] We note that partial wave unitarity of e~e~ — > H~ H~ also sets 
a bound on £ which is, however, weaker than £xj coming from 
Tat < mjv/2. We have also checked that the measured g — 2 of 
the electron does not further constrain £ . 



GeV ee machines. For illustration, we take to#+ = 150 
and 300 GeV, for the ^ = 500 and 1000 GeV colliders, 
respectively. We consider the cases £ = 0.1 and £ = 1, 
imposing the unitarity bound such that £ = min(0.1, 
and £ = min(l, £ry), respectively (recall that the unitarity 
bound implies £ry ~ 0.1 already for tojv ~ 7 TeV). As 
expected, there is a sharp decrease in a which occurs 
once £ enters the unitarity bound regime (£ — > £{/), for 
which the cross-section drops as to^ 2 as seen in ([9]). 

In the mass range 100 GeV & m N < few • 10 3 TeV, de- 
spite fixing £ = min(0.1,£ry) due to unitarity, we still see 
that a ~ 10~ 3 • oqed- Thus, recalling that any given 
design of an ee machine relies on having 10000 events per 
unit (Jqed, one expects more than 10 e~e~ — > H~H~ 
events even if to^ ~ 10 3 — 10 4 TeV. In Fig. [2] we plot 
a "naive" discovery limit (demanding 10 events for dis- 
covery) for e~e~ — > H~H~ in the — £ 2 plane, at 
an ELC with y/s = 500 GeV, 1 TeV and 3 TeV, setting 
m,H = 150, 300 and 600 GeV, respectively, and with an 
integrated luminosity of L = lOOx (y/s/TcV) 2 fb _1 which 
corresponds to L ■ <jqed ~ 10000 at each cm. energy. 
We see that the like-sign Higgs pair production signal is 
indeed striking at an ILC; it is potentially observed even 
if rriM ~ 10 4 TeV and £ 2 ~ 10 -8 - an extremely small 
value for the effective neutrino mixing angle which is con- 
sistent with the unitarity bound and more than 7 orders 
of magnitudes smaller than the (3(3q v bound at to at ~ 10 4 
TeV! 




m N [GeV] 



FIG. 2: Discovery limit for e~ e~ — > H~H~ (long-dashed 
lines) and for e~e + — ► H~ H + (dotted lines), at an ILC with 
y/s — 500 GeV (upper line), y/s — 1 TeV (middle line) and 
y/s = 3 TeV (lower line), see also text. The parameter space 
above the lines corresponds to the observable events. Also plot- 
ted (solid lines): the experimental limits on U^ N from (3(3q v 
and on |£/ e jv| 2 from invisible Z-decays, as well as the perturba- 
tive unitarity bound from Q). Here the parameter space above 
the lines are ruled out. 

Let us now turn to the lepton number conserving 
(LNC) process e + e~ — > H + H~ . There are 3 diagrams 
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contributing to this process: the two "standard" (i.e., 
within multi-Higgs models) s-channel 7 and Z exchanges 
and a 3rd t-channel N exchange one. The cross-section 
a h = a[eT ~e + — > H~H + ) can be decomposed accord- 
ingly as a h = a~ s + + a^ + cr~y , where the subscripts 
S, N and SN denote the standard s-channel 7, Z con- 
tribution, the t-channel N contribution and the inter- 
ference between the standard and the V-exchangc dia- 
grams. While as can be found in the literature [18|], the 
new N exchange contributions to the total cross-section 
are: 



a SN = °n 



1 - 4r z s w c w 



'N 



= c 



N- 



2e 



P 2 r w 



In 
w In 



w — (3 
w + /3 



■0 



f3w 
2/3 



where v = rjv + (1 — w) 2 /4 and 



Cn = cfqed ■ £ 



2 3v rN 
32s^ r w 



,(10) 



(11) 
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FIG. 3: The cross-sections ratios a /c(e e — > W W 
and a l ~/cr(e _ e + — > W~W + ), as a function of uln, for £ 
UeN = min(l, £u) at a 500 GeV ILC. 



For a constant £ and in the limit m 2 N >> s, both 
tJg^ and a^ + approach a constant with respect to m^. 
As £ enters the perturbative unitarity bound regime, 
i.e., £ — > £[/ ~ 1/tojv, the interference term dominates 
SN/ a N + m Af/ s , b u t still drops with tojv as 



since (7 



^SJV ^ l/ m Af- Thus, in this large N mass range the to- 



cross-section converges to the standard value: 



taliT 1 

a l ~(mf v >> s,£ = £u) — * a s + ■ This behavior can be 
seen in Fig. [TJ where we have plotted a h in units of 
the point cross-section uqed, as a function of for 
both a 500 and 1000 GeV ee machines, for which we take 
m H + — 150 and 300 GeV, respectively. As in the LNV 
H~H~ case, we consider £ = 0.1 and £ = 1, impos- 
ing the unitarity bound such that £ = min(0.1,£j/) and 

£ = nun(l,£tf)- 

Note that around ~ ^/i the t-channel N exchange 
in H + H~ production can become sizable. We, thus, de- 
fine the number of standard deviations (Ng D ) with which 
the t-channel N exchange in e~e + — > H~H + can be de- 
tected at an ILC: 



N N - 

Iy SD — 



V<F+ 



(12) 



In Fig. [2] we plot the 5a (Ng D — 5) discovery contours 
of the V-exchange contribution to e~e + — > H~H + in the 
£ 2 — ton plane. As in the H~H~ case, we take an ILC 
with V^s = 500 GeV, 1 TeV and 3 TeV (setting m H = 
150, 300 and 600 GeV, respectively), and an integrated 
luminosity of L = 100 x (^/i/TeV) 2 fb _1 . Although not 
as promising as the LNV H~H~ signal, we see that the 
N- mediated H~H + LNC signal can be detected at 5a 
in an ILC with e.g., s/s = 1(3) TeV, for a very heavy N 
with a mass in the range 10(2) TeV < m N < 30(100) TeV 
and £ 2 ~ 6 • 10~ 4 (6 • 10~ 5 ) - a value well within the (3f3 0u 
bound and consistent with the unitarity bound. 



Admittedly, our criteria for a 5a discovery (or 10 events 
in the H~H~ case), as plotted in Fig. [21 is rather naive 
since we do not take into account the possible charged 
Higgs decays, their branching ratios and the correspond- 
ing detection efficiencies. Nonetheless, this simplified dis- 
covery criteria allows one to get an idea of the rather 
striking prospects of these charged Higgs pair produc- 
tion signals (especially the LNV H~H~ one), as po- 
tential probes of heavy Majorana neutrinos. In par- 
ticular, to further appreciate the significance of both 
the LNV and LNC Higgs pair production signals, we 
plot in Fig. [3] the ratios a /cr(e~e~ — > W~W~) and 
a-+/a(e~e+ -► W~W + ) at a 500 GeV iLC,^ 3 ] as a 
function of mjv and for £ = U c n = mm(l, £(/), i.e., as- 
suming = 1. 

Take for example the H~H~ signal for which we get 
a ja(ere~ — > W~W~) > 10 already for tun }t y/s 



(=500 GeV). Recall that a(e 



W~W~) oc U\ 



N 



the inverse /30ov (t-channel N exchange) process, which 
was extensively studied as a potential discovery chan- 
nel for a heavy N at an ILC [12], and which stands 
as the leading background to the e~e~ — * H~H~ pro- 
cess. Thus, the H~H~ signal should be easily de- 
tected at a 500 GeV ILC if m N £ 500 GeV. Clearly, if 
mjj > m t + mb, then the H~ H~ signal can be searched 
for via e - e~ — > H~H~ — > bbtt. In this case the effective 
H~ H~ /W~W~ ratio is further (dramatically) enhanced 
by a (2 — > 2)/(2 — > 4) phase-space factor. As for the LNC 
H-H+ channel, we get a~ + /a{e-e + -> W~W + ) > 0.1 
for 100 GeV < m N < 1000 GeV. Thus, here also, in par- 
ticular if run > nrit the N-exchange contribution in 
the H~H + channel should be easily accessible to a 500 
GeV ILC if m N - 0(1) TeV, even if £ 2 < 10" 2 - see also 

Fig. m 

To summarize, we have examined the potential role 
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of the charged Higgs pair production processes e~e + — > 
H~H + and e~ e~ — > H~H~ at an ILC, as possible 
probes of new dynamics mediated by heavy neutrinos, 
N. We have shown that the mjsi/mw enhancement fac- 
tor in the e — N — H + interaction vertex (i.e., compared 
to the e ~ N — W + vertex) makes these processes ex- 
tremely sensitive to the t-channel iV-exchange that drives 
them. In particular, we find that the LNV H~H~ sig- 
nal is sensitive to a heavy neutrino with a mass up to 
10 4 TeV and an effective Yukawa-like coupling to an 
electron as small as £ ~ 1CP 4 . For example, in a su- 
persymmetric non-minimal seesaw-like framework with 
heavy Major ana neutrinos, for which £ = U & n ■ and 
= 1/ tan/3, one expects £ = U e N / tan/3, where U & n is 
the v e — N mixing angle. Thus, taking tan/3 ~ O(30) or 
tan/3 ~ 0(3), the LNV H~H~ signal can probe a very 
heavy N with win ~ 10 4 TeV for the very small values of 
the neutrino mixing angles U^ N ~ 1CT 5 or U* N ~ 10 ~ 7 , 
respectively. This striking result is compatible with both 
the /3/3o„ and the perturbative unitarity bounds on the 
UeN~ m N plane and, to the best of our knowledge, stands 
out as an exceptional signal of heavy Majorana neutrinos 
at large colliders. 

In the case of the LNC signal e + e~ — ► H + H~ , we find 
that an ILC with a center of mass energy of 1 TeV, can 
potentially probe a heavy Majorana with a mass mjv ~ 
O(10) TeV, if £ 2 - 0.001 - also within the /3/3 „ and the 

[2] We have included the contribution from the t-channel N- 

exchange diagram in a(e~e + — > 
[3] Following suggestions by Tao Han, we have also considered the 

constraints on t/ e jv from the LEP2 measurements of a(e~e + — > 



perturbative unitarity bounds on the U^ N — win plane. 

Finally, in passing, we note that recently an interesting 
model that can be tested at the ILC through our charged 
Higgs pair production mechanisms was suggested by Ma 
in [3J. In Ma's model a second scalar doublet (r) + ,rf) is 
introduced and the light neutrinos acquire sub-eV masses 
through a loop involving both N and rf - a radiative 
seesaw mechanism. As a consequance, the coupling e — 
N — rf~ is much less constrained (than in the minimal 
seesaw scenario) and should, therefore, be studied at the 
ILC via the reaction e~e~ — ► r\~rf~ , where each of the 
?/ _ 's in the final state can further decay via rf~ — > rfW~, 
rf being the SU(2) partner of r/ + . It is also worth noting 
that, in this model, rf can be the lightest stable particle 
due to an exact Zi symmetry of the model, making it a 
good candidate for dark matter ||. In such a case, the 
reaction e - e~ — ► T]~r}~ will lead to a unique signal of 
two same-charge IF-bosons and a missing energy carried 
by the two rj°'s. 

We thank Gad Eilam, Tao Han, Ernest Ma and Jose 
Wudka for very useful discussions. S.B.S thanks the hos- 
pitality of the theory group in Brookhaven National Lab- 
oratory where part of this study was performed. This 
work was supported in part by US DOE Contract Nos. 
DE-FG02-94ER40817 (ISU) and DE-AC02-98CH10886 
(BNL). 

W~W+) 0. We found that e"e+ -» W~W+ is not sensitive 
at LEP2 to the N-exchange diagram if U e M 0.1, i.e., within its 
bound from invisible Z decays. 
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